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ABSTRACT 

The  US  Naval  Research  Laboratory  (NRL)  is  developing  technolo¬ 
gies  that  will  enable  Navy-relevant  missions  with  the  smallest  practi¬ 
cal  Micro  Air  Vehicles  (MAVs).  The  NRL  Micro  Tactical  Expend¬ 
able  (MITE)  air  vehicle  is  a  result  of  this  research.  MITE  is  a 
hand-launched,  dual-propeller,  fixed-wing  air  vehicle,  with  a  25cm 
chord  and  a  wingspan  of  25 -47cm,  depending  on  payload  weight. 
Vehicle  gross  weight  is  130-350g.  Miniature  autopilot  systems, 
based  on  visual  imaging  techniques,  are  being  developed  for  MITE. 
These  will  be  used  in  conjunction  with  conventional  autopilot  sen¬ 
sors  to  allow  the  MITE  to  fly  autonomously.  This  paper  provides  an 
overview  of  the  MITE  development,  including  aerodynamic  design 
considerations,  electric  propulsion,  and  vision-based  autopilot  re¬ 
search.  Also  presented  is  a  rationale  for  the  development  of  control 
laws  that  can  direct  the  behavior  of  large  groups  of  MAVs  or  other 
vehicle  agents.  Dubbed  ‘physicomimetics,’  this  process  can  bring 
about  the  self-assembly  of  complex  MAV  formations,  though  indi¬ 
vidual  MAVs  have  minimal  onboard  processing  power  and  limited 
local  sensing  capabilities. 


NOMENCLATURE 

AR  aspect  ratio 

b  wing  span 

c  wing  chord 

Cd  total  drag  coefficient 

Cdi  induced  drag  coefficient 

Cd  0  profile  drag  coefficient 

Q  lift  coefficient 

Cmx  roll  moment  coefficient 

Cmz  pitching  moment  coefficient 

e  lift  distribution  efficiency 

F  force 

G  physicomimetic  gravitational  constant 


Kz  constant 

L  characteristic  length  (wing  chord) 
M  particle  mass 
p  particle  position,  (x,y) 
r  range 
R  radius 

Re  reynolds  number,  t/Lp/p 
t  time 

a  angle-of-attack 
\x  air  viscosity 
v  velocity 

0  differential  elevon  deflection  angle 
p  air  density 


INTRODUCTION 

Micro  Air  Vehicles  (MAVs)  offer  the  promise  of  affordably  expend¬ 
able,  highly  portable  covert  sensor  platforms  for  a  range  of  close-in 
situational  awareness  activities.  Since  1996,  the  US  Naval  Research 
Laboratory  (NRL)  has  been  sponsored  by  the  Office  of  Naval  Re¬ 
search  to  develop  technologies  that  will  enable  Navy-relevant  appli¬ 
cations  with  the  smallest  practical  MAVs.  The  goal  of  this  research 
is  to  develop  and  demonstrate  MAV  technology  that  supports  Navy 
and  Marine  Corps  missions  such  as  over-the-hill  reconnaissance, 
area  or  perimeter  patrol,  and  battle  damage  assessment.  NRL’s  work 
is  both  complimentary  and  supplementary  to  the  Defense  Advanced 
Research  Projects  Administration  (DARPA)  MAV  Program.  Addi¬ 
tional  MAV  requirements  beyond  the  DARPA  effort  include  the 
Navy-specific  applications,  electric  propulsion,  non-Global  Position¬ 
ing  System  (GPS)  navigation,  and  airframe  size  commensurate  with 
operating  conditions. 
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Figure  1 . 


MITE  VEHICLE  DESIGN 

NRL  developed  a  baseline  design  called  the  Micro  Tactical  Expend¬ 
able  (MITE)  which  consists  of  a  low  aspect  ratio  flying  wing  with 
dual,  counter-rotating  propellers  mounted  at  the  wing  tips  (see  Fig. 
1).  This  configuration  was  selected  based  on  a  tradeoff  analysis  that 
included  aerodynamic  performance,  payload  interface,  launch  and 
recovery  methods,  and  compact  storage.  Specifically,  the  low  wing 
aspect  ratio  is  necessary  to  provide  sufficient  wing  area  within  a 
compact  wingspan.  Although  the  low  aspect  ratio  increases  the  in¬ 
duced  drag  of  the  wing,  the  increase  in  wing  chord  raises  the 
Reynolds  number,  which  improves  the  boundary  layer  characteris¬ 
tics,  hence  performance,  of  the  aerofoil.  The  dual  propellers  provide 
slipstream  flow  over  nearly  the  entire  wing  for  enhanced  lift  at  low 
speeds.  By  counter  rotating  in  the  direction  opposite  of  the  wingtip 
vorticies,  the  Zimmerman  Effect^)  is  produced,  which  reduces  the 
induced  drag  to  a  value  noticeably  lower  than  expected  for  the  low 
aspect  ratio  of  the  wing.  Another  benefit  of  propeller  counter-rota¬ 
tion  is  the  balancing  of  torque  and  slipstream  rotation  effects,  allow 
ing  high  lateral  stability  at  low  speeds  to  enable  easy  hand  launching 
of  the  MITE.  As  a  flying  wing,  MITE  uses  elevon  controls,  which 
have  enhanced  effectiveness  since  they  are  also  immersed  in  the  slip¬ 
stream.  Payload  location  is  in  an  unobstructed  nose,  which  is  ideal 
for  imaging  and  accommodating  a  range  of  sensors.  Despite  being  a 
flying  wing  with  no  geometric  dihedral,  the  overall  combination  of 
design  features  results  in  a  simple  configuration  possessing  inherent 
stability  about  all  three  axes,  good  performance,  ease  of  hand  launch, 
and  simple  landing  by  gliding  to  the  ground. 

Wing  aspect  ratio  at  low  Reynolds  numbers 

Fifty  percent  of  the  total  drag  for  a  well-designed  aeroplane,  cruising 
at  its  best  speed  for  range,  is  induced  drag.  When  flying  at  best  effi¬ 
ciency  for  endurance,  the  induced  drag  component  is  an  even  greater 
percentage  of  the  total  drag.  The  profile  drag  of  the  aerofoil  shape  is 
a  relatively  small  component  of  the  overall  aeroplane’s  drag.  Conse¬ 
quently,  while  aerofoil  designers  strive  to  minimise  profile  drag, 
aeroplane  designers  must  minimise  induced  drag  to  maximise  range 
and  endurance. 

The  induced  drag  of  a  wing  is  a  result  of  a  spanwise  component  of 
the  airflow  near  the  wingtips  caused  by  the  finite  span  of  a  wing.  In¬ 
duced  drag  is  a  function  of  the  square  of  the  lift  coefficient  and  in¬ 
versely  proportional  to  the  wing  aspect  ratio.  Since  the  optimal 
cruise  lift  coefficient  is  a  function  of  the  aerofoil,  reducing  induced 
drag  requires  increasing  the  aspect  ratio,  i.e.  increasing  the  wing 
span-to-wing  chord  ratio.  For  a  given  wing  loading,  increasing  the 
wing  aspect  ratio  results  in  a  reduction  of  the  wing  chord.  Typically, 
the  wing  aspect  ratio  is  as  large  as  structural  considerations  allow, 
since  reducing  the  chord  reduces  the  wing  thickness  in  proportion, 


Span 

Configuration  A 

36-8cm 

Configuration  B 

36-8cm 

Gross  Weight 

129  grams 

211  grams 

Motors 

two  4  watt  coreless,  geared 

two  7  watt 
coreless,  geared 

Porpellers 

17- 8cm  dia,  20- 3cm  pitch 

17-8cm  dia,  20-3cm 
pitch 

Propulsion  Battery 
(CR2) 

9  volt  LiSQ  primary  (CR2) 

12  volt  LiSQprimary 

Command  Receiver 

single  conv.  FM  72MHz 

dual  conv  FM 

72MHz 

Payload 

none 

colour  video  system 

Table  1 

while  at  the  same  time  increasing  the  wing  root  bending  moment 
since  the  wingspan  is  increased. 

Unlike  large  aircraft,  the  maximum  value  of  wing  aspect  ratio  for 
an  MAV  is  limited  by  aerofoil  aerodynamics  before  the  structural 
limit  is  reached.  The  value  of  profile  drag  for  a  large  aircraft  is  rela¬ 
tively  constant  with  Reynolds  number  variations.  However,  in  the 
aerodynamic  regime  of  MAV  aerofoils,  the  profile  drag  coefficient 
increases  in  inverse  proportion  to  the  cube  root  of  the  Reynolds 
number.  As  the  value  of  the  aerofoil’ s  profile  drag  begins  to  rise,  a 
point  is  reached  where  it  is  no  longer  a  minor  drag  component.  As 
such,  the  design  of  an  efficient  MAV  wing  requires  an  optimisation 
between  minimising  induced  drag  without  significantly  increasing 
the  aerofoil  profile  drag.  An  empirical  equation  (Fig.  2)  was  derived 
at  NRL  that  can  be  used  to  determine  an  optimal  tradeoff 
between  induced  drag  and  profile  drag  for  MAVs.  This  Equation 
substantiates  why  MAVs  typically  perform  well  with  wings  having 
lower  aspect  ratio  than  larger  aircraft.  For  MITE,  chord  Reynolds 
number  ranges  between  80,000  and  200,000. 

MITE  prototype  development 

A  number  of  exploratory  MITE  flight  test  vehicles  have  been 
designed  and  built  to  explore  various  vehicle  design  configurations, 
including  a  range  of  propulsion  and  control  approaches.  Design 
focused  on  minimum  size  air  vehicles  that  were  capable  of  perform 
ing  potentially  useful  military  missions.  Maximum  use  was  to  be 
made  of  available,  commercial  items,  while  exploiting  emerging, 
high  pay-off  technology  where  appropriate.  The  vehicles  were  also 
available  to  support  the  initial  flight-testing  of  prototype  and  bread¬ 
board  payload  components  and  modules.  The  test  vehicles  ranged  in 
size  from  25 -47cm,  with  gross  weights  of  130-350g. 

Early  tests  of  miniature  internal  combustion,  liquid  fuel  engine 
propulsion  were  abandoned  in  order  to  concentrate  on  electric 
motors  employing  high  energy-density  primary  (single  use)  and  sec¬ 
ondary  (rechargeable)  cells.  Electric  propulsion,  while  not  as 
weight/energy  efficient  as  chemical  fueled  internal  combustion 
engines,  offered  the  potential  for  highly  reliable,  totally  silent  field 
operations  plus  extremely  rapid  commercial  technological  growth. 
Early  flight-testing  of  various  motor  configurations  also  indicated 
that  twin  motor  layouts,  compared  to  single  motors,  could  provide 
simplified  launching  capability  and  possible  control  and  aerodynam 
ic  advantages.  Plank- type  flying  wing  configurations  proved  to  be 
compact,  sturdy,  and  stable  MITE  platforms.  Elevon  plus  synchro¬ 
nised  or  differential  motor  speed  control  provides  adequate  vehicle 
control.  The  principal  MITE  vehicles  are  described  below. 

MITE  2  (see  Fig.  3),  the  first  twin  motor  flight  test  vehicle,  was 
flown  in  two  configurations,  as  seen  in  Table  1.  The  MITE  2  air- 
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Cdo  =  Ki 

1/3 

(Re) 


Where  KI  is  determined 
from  empirical  airfoil  data 

Example:  For  the  Selig  SI  091 
Cdo  =  0.0165(60,000/Re)1/3 


Cd  —  Cdo  +  Cdi 
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The  Profile  Drag  Coefficient  For  Low  Reynolds  Number 
Airfoils  Varies  Inversely  Proportionally  To  The  Cube 
Root  Of  Reynolds  Number 


The  Wing's  Drag  Is  The  Sum  Of  The  Profile  Drag  Of  Th\ 
Airfoil  And  The  Induced  Drag  Due  To  Having  A  Finite 
Wing  Span 


Chord  Reynolds  Number  Equals  The  Product  Of  The 
Air  Density,  The  Airspeed,  And  The  Chord  Length 
Divided  By  The  Viscosity 

Aspect  Ratio  Is  The  Wing  Span  Divided  By  The 
Average  Chord 


So, 


^  Ki  /  K3  K2  /  b 
(C)1/3  (1/c) 

=  Ki'(c)1/3  +  K2'(C) 

dCD/  dc  =  0  =  -1/3  Ki'(c)4/3+  K2' 


Where  Ki'  =  Ki/K3 
& 

K2'  =  K2  /  b 


For  Minimum  Drag  dCo/  dc  =  0 


3/4 

3/neb 

Wing  Chord  For  Minimum  Wing  Drag 

Therefore, 

c  = 

.Ki  C^c//ov)1/3_ 

Figure  2. 


frame  structure  utilised  balsa,  foam,  carbon  epoxy,  and  a  polyester 
fiber  covering.  In  early  tests  using  lithium  primary  cells,  it  demon¬ 
strated  a  maximum  endurance  of  about  30  minutes.  It  was  a  stable, 
and  relatively  smooth  camera  platform.  A  competent  pilot  could  fly 
the  vehicle  via  the  on  board  video  system,  operating  at  airspeeds  of 
4-8  ms- !  Valuable  experience  was  gleaned  from  numerous  flights 
with  an  on-board  video  system.  Typical  weight  distributions  of 
MITE  2  vehicles  are  shown  in  Table  2. 

MITE  3  (see  Fig.  4)  was  a  30cm  span  version  of  the  MITE  2B.  It 
incorporated  more  powerful  control  system  servos  and  electronic 
speed  controls  (ESC).  The  MITE  3  and  subsequent  versions  also  em 
ployed  lithium  ion  secondary  cells  for  propulsion  batteries  because 
of  their  superior  energy  and  power  density  at  relevant  power  levels. 
These  batteries  provided  up  to  9  watts  output  per  motor.  Gross 
weight  without  payload  was  129g. 

MITE  4  (see  Fig.  5),  the  most  recent  configuration,  is  a  47cm  span 
vehicle  designed  to  be  able  to  carry  a  variety  of  early  prototype  or 
breadboard  payload  components  of  85-1 15g  for  test  purposes.  In  ad¬ 
dition  to  the  larger  wing  area  (l,200cnf  vs  650cm?  for  MITE  3),  it 
carries  somewhat  more  powerful  (approximately  10  watts)  geared 
coreless  motors.  Gross  weight  without  payload  is  about  210-260g, 
depending  on  the  propulsion  battery  pack.  Payloads  of  up  to  lOOg 
can  be  accommodated 

Motor/battery  system 

The  electric  motors  used  in  ah  MITE  vehicles  to  date  are  commer¬ 
cially  available,  coreless,  permanent  magnet  motors.  At  relatively 
low  output  power  levels  of  under  15  watts,  coreless  motors  are 
lighter  in  weight  and  generally  more  efficient  than  conventional  fer¬ 
rite  cored  motors.  Their  electrical  efficiency  can  be  above  70%,  al¬ 


though  at  power  levels  usually  employed,  60%  efficiencies  are  more 
common.  Because  these  motors  are  high  RPM,  low  torque  devices, 
gearing  is  used  in  ah  cases  to  turn  efficient  size  propellers  at  lower 
speeds  than  the  motor  rotor.  Gearing  of  up  to  12:1  has  been  used,  but 
current  MITE  vehicles  typically  use  around  6:1  gear  ratios  to  turn 
17-8cm  diameter,  20-3cm  pitch  propellers  at  speeds  between  3,000 
and  5,000rpm.  To  achieve  symmetrical  propeller  flow  forces  on  the 
vehicle,  opposite  rotation  is  used  on  the  two  motors.  This  allows  eas¬ 
ier  launching,  including  hand  launching,  of  the  vehicles  compared  to 
single  propeller  configurations.  By  turning  the  propellers  such  that 
the  inner  portions  of  the  blades  (closest  to  the  fuselage)  are  moving 
upward,  it  is  also  possible  to  gain  some  aerodynamic  benefits.  Bat¬ 
teries  tested  for  use  in  the  MITE  vehicles  included  nickel  cadmium 
(NiCd),  nickel  metal  hydride  (NiMH),  and  lithium  ion  secondary 
cells  as  well  as  lithium  sulfur  dioxide  (LiS02)  primary  cells.  For  dis¬ 
charge  rates  of  1  ampere  or  less,  the  lithium  primary  cells  proved  to 
be  superior  to  all  other  cells  tested.  The  nominal  3 volt,  750mAh  ca¬ 
pacity  CR2  cell,  weighing  less  than  12g,  was  most  efficient  for  use 
in  the  early  lightweight  MITE  2  test  vehicles.  For  subsequent  MITE 
vehicles,  where  higher  power  and  current  drains  of  3-4  amperes  was 
needed,  lithium  ion  secondary  cells  of  430  and  800mAh  nominal  ca¬ 
pacity  were  used.  As  shown  in  Table  3,  these  cells  provide  much 
higher  energy  densities  than  comparable  weight  nickel  cadmium  or 
nickel  metal  hydride  cells. 

Design  refinement  via  computational  fluid  dynamics 

Computational  fluid  dynamics  (CFD)  plays  an  important  role  in  the 
design  of  low  Reynolds  number  UAVs  or  MAYs.  NRL  has  been  us¬ 
ing  the  Finite  Element  Flow  solver  code  (FEFLOy2),  an  unstructured 
finite  element  code  based  on  tetrahedral  elements,  to  investigate  a 
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Figure  3. 
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Figure  4. 


range  of  small  vehicles  suitable  for  autonomous  flight  operation. 
CFD  is  used  for  two  objectives.  One  is  to  determine  optimal  config¬ 
urations  from  the  aerodynamic  performance  point  of  view.  The  other 
is  to  determine  the  aerodynamic  coefficients  that  serve  as  input  to 
NRL’s  six  degrees  of  freedom  flight  simulator.  This  simulator  is  in 
turn  used  to  evaluate  stability,  performance,  and  optimal  control 
laws  for  the  autopilot. 

The  FEFLO  code  was  developed  at  NRL  by  Ravi  Ramamurti  and 
at  The  George  Mason  University,  Fairfax,  Virginia,  by  Rainald 
Lohner.  The  governing  equations  are  the  incompressible  Navier- 
Stokes  Equations,  which  are  solved  in  an  Arbitrary  Lagrangian 
Eulerian  (ALE)  formulation.  The  Equations  are  discretised  in  time 
using  an  implicit  time  stepping  procedure.  It  is  important  for  the 
flow  solver  to  be  able  to  capture  the  unsteadiness  in  the  flow  field,  if 
such  exists.  This  is  especially  true  in  low  Reynolds  number  vehicles 
where  the  separation  bubble  may  be  unsteady.  The  flow  solver  is 
time- accurate,  allowing  local  timestepping  as  an  option.  The  result¬ 
ing  expressions  are  subsequently  discretised  in  space  using  a 
Galerkin  procedure  with  linear  tetrahedral  elements.  In  order  for  the 
algorithm  to  be  as  fast  as  possible,  the  overhead  in  building  element 
matrices,  residual  vectors,  etc.  must  be  kept  to  a  minimum.  This  is 
accomplished  by  employing  simple,  low-order  elements  that  have  all 
the  variables  (three  velocity  components  and  pressure)  at  the  same 
node  location.  The  resulting  matrix  systems  are  solved  iteratively  us¬ 
ing  a  preconditioned  conjugate  gradient  algorithm  (PCG).  The  code 
has  the  options  to  be  run  Euler,  laminar  flow,  or  incorporate  a  Bald- 
win-Lomax  turbulence  model.  The  solver  has  been  successfully 
evaluated  for  both  two  dimensional  (2D)  and  three  dimensional 
(3D),  laminar  and  turbulent  flow  problems  by  Ramamuti  et  al (3  and4). 

Flow  over  small  vehicles  is  significantly  modified  by  the  influ¬ 


ence  of  the  propeller.  Time  dependent  simulation  of  the  propeller 
would  be  very  time  consuming,  and  inconsistent  with  the  desire  to 
make  many  simulations  of  different  configurations,  at  several  vehi¬ 
cle  attitudes,  to  determine  the  aerodynamic  coefficients.  Instead,  it 
has  been  found  to  be  adequate  to  model  the  propeller  as  an  actuator 
disk  with  a  distribution  of  axial  and  radial  momentum  sources  that 
match  the  thrust  and  drag  of  the  actual  propeller. 

Figure  6  shows  the  flow  field  over  two  different  configurations. 
The  first  (Fig.  6(a))  is  the  MITE  vehicle  with  a  wingspan  of  15cm 
and  a  span  to  wing  cord  aspect  ratio  of  1-25.  The  pressure  contours 
on  the  body  surface  are  shown.  The  dark  regions  indicate  stagnation 
regions  of  high  pressure.  The  high  pressure  in  the  stagnation  regions 
around  the  wing  body  junction  can  be  reduced  by  proper  fairing.  A 
number  of  configuration  changes  were  looked  at  to  find  an  optimal 
configuration.  The  inviscid  calculations  were  performed  with 
approximately  160k  points  and  860k  tetraheda.  For  symmetric  cases 
the  simulations  were  performed  over  half  of  the  vehicle  assuming  a 
symmetry  plane  at  the  axis.  Figure  6(b)  shows  the  MITE  2  configu¬ 
ration,  a  slightly  larger  vehicle  with  a  wingspan  of  36-8cm  and  a 
cord  of  254cm.  Also  shown  are  the  velocity  contours  for  a  simula¬ 
tion  that  included  effects  of  the  counter  rotating  propellers.  The  pro¬ 
peller  rotation  direction  was  chosen  to  counter  the  effect  of  the  tip 
vortex  and  increase  the  effective  aspect  ratio.  This  simulation  was  at 
an  angl-  of-attack  of  15°.  Note  the  induced  separation  at  the  wing 
body  junction  on  the  MITE  2  configuration  at  the  high  angle  of  at¬ 
tack.  The  results  of  the  aerodynamic  coefficients  derived  for  the 
MITE  2  configurations  are  shown  in  Fig.  7. 


Table  2. 

MITE  2  Weight  Breakdown 


Item  Weight,  grams 


Ver  2A 

Ver  2B 

Command  receiver 

5 

21 

Control  servos 

6 

6 

Motor  speed  control 

5 

6 

Motors 

25 

25 

Props 

7 

7 

Battery  power  supply 

34 

45 

Cable  Harness 

6 

8 

Airframe  structure 

41 

41 

Sub-total 

129 

159 

Payload 

0 

52 

gross  weight 

129 

211 

Figure  5. 
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Cell 

Weight, 

Discharge 

Energy 

Energy 

AV 

Delivered 

Cutoff 

Cutoff 

g 

rate, 

capacity 

density, 

voltage 

capacity 

voltage 

time,  min 

amps 

watt  min 

watt 

min/g 

mAh 

430  mAh  Tadiran 

12-5 

0-7 

59-2 

4-73 

2-62 

376 

2-2 

32-2 

Lithium  Rechargeable 

TO 

44-9 

3-59 

2-53 

294 

2-2 

17-7 

T5 

32-0 

2-56 

247 

216 

2-2 

8-6 

2-0 

18-3 

147 

2-36 

130 

2-2 

3-9 

800  mAh  Tadiran 

18-3 

TO 

93-0 

5-08 

2-55 

608 

2-2 

36-5 

Lithium  Rechargeable 

T5 

82-6 

4-50 

2-56 

538 

2-2 

2T5 

2-0 

7T7 

3-92 

2-51 

476 

2-2 

14-3 

2-5 

60-2 

3-29 

2-39 

420 

2-2 

104 

3-0 

46-7 

2-55 

240 

324 

2-2 

6-5 

750  mAh  CR2 

11-6 

0-7 

82-9 

745 

2-30 

595 

T6 

5T0 

Lithium  Pimary 

TO 

57-1 

4-92 

2-20 

425 

T6 

25-0 

T5 

34-6 

2-98 

240 

280 

1.6 

10-3 

50  mAh  Sanyo  NiCd 

3-6 

0-5 

4-7 

T31 

144 

35 

0-8 

44 

TO 

T8 

0-50 

0-99 

31 

0-8 

T8 

T5 

T1 

0-29 

0-99 

27 

0-8 

14 

80  mAh  Sanyo  NiCd 

5-73 

0-5 

8-7 

T53 

147 

63 

0-8 

7-5 

TO 

3-3 

0-57 

T09 

50 

0-8 

3-0 

T5 

T8 

0-32 

T03 

44 

0-8 

T8 

110  mAh  Sanyo  NiCd 

7-3 

0-5 

6-7 

0-91 

149 

94 

0-8 

1T3 

0-7 

6-7 

0-91 

146 

90 

0-8 

7-8 

TO 

6-1 

0-82 

142 

90 

0-8 

54 

T5 

5-5 

0-74 

T07 

86 

0-8 

3-5 

2-0 

5-2 

0-70 

TOO 

82 

0-8 

2-5 

120  mAh  NiMH 

3-5 

0-7 

7-1 

2-00 

T01 

105 

0-8 

94 

TO 

5-7 

T60 

0-95 

92 

0-8 

5-6 

T5 

4-9 

140 

0-94 

84 

0-8 

34 

270  mAh  NiMH 

7-5 

T5 

17-0 

2-30 

T03 

278 

0-8 

114 

2-0 

15-0 

2-00 

0-95 

255 

0-8 

7-7 

2-5 

13-0 

T70 

0-93 

234 

0-8 

5-6 

3-0 

10-0 

T30 

88 

135 

0-8 

3-7 

Table  3. 

Sample  Cell  Test  Comparison 


MITE  AUTOPILOT  DEVELOPMENT 

Autopilots  provide  a  significant  challenge  for  small  autonomous 
vehicles.  While  there  are  many  UAV  autopilots  in  production  for 
operational  systems,  most  commercial-off-the-shelf  autopilots  weigh 
2-6  Kg.  This  is  much  greater  than  the  gross  weight  of  the  MITE 
class  of  vehicles.  NRL  has  developed  a  number  of  autopilots  in  sup¬ 
port  of  small  UAVs.  A  recent  example  is  the  Lightweight  Autopilot 
(LWAP),  which,  at  900g  meets  the  functional  but  not  the  physical 
requirements.  In  an  attempt  to  reduce  the  weight  of  the  autopilot  we 
are  pursuing  a  number  of  unconventional  means  of  autopilot  control. 
These  are  primarily  aimed  at  reducing  the  number  of  sensors 
required  to  provide  the  necessary  stability  and  directional  control.  To 
test  these  new  algorithms,  a  six  degree  of  freedom  (6  DOE)  flight 
simulator  was  used,  coupled  to  an  autopilot  routine  that  uses  simulat¬ 
ed  sensors  with  realistic  resolution  and  noise.  The  CLD  simulations 
are  used  to  provide  the  aerodynamic  coefficients  and  stability  deriva¬ 
tives  for  the  6  DOE  code.  Two  types  of  algorithms  are  being 
pursued.  One  involves  using  data  from  optical  sensors  to  provide 
attitude  information  or  optical  flow  for  collision  avoidance.  The 
other  uses  a  more  sophisticated  processing  of  a  more  limited  set  of 
sensor  input  and  relies  on  the  inherent  stability  of  the  MITE  2  vehi¬ 
cle  to  achieve  directional  and  altitude  control.  A  control  law  has 
been  found  that  will  allow  the  MITE  2  vehicle  to  be  flown 
autonomously  with  only  altitude  and  heading  information  that  can  be 
provided  by,  for  instance,  a  simple  barometric  pressure  sensor  and 
magnetic  heading  indicator.  Additional  sensors  can  be  added  that 
increase  the  functionality,  albeit  at  the  expense  of  additional  weight. 

Optical  flow  sensing  for  collision  avoidance 


Optical  flow  sensing  is  a  technique  that  allows  a  moving  observer  to 
sense  the  proximity  of  its  surroundings  and  relative  motion  with  a 
minimum  of  processing  power.  Flying  insects,  though  equipped  with 
only  a  rudimentary  vision  system,  efficiently  use  optical  flow  to  es¬ 
cape  attack,  avoid  obstacles,  and  to  alight  on  objects.  An  optical 
flow  sensor  does  not  produce  an  image  as  a  camera  or  a  human  eye 
would.  Rather,  it  senses  the  relative  movement  of  the  texture  of  its 
surroundings  by  tracking  contrasting  regions  of  light  and  shade  as 
the  sensor  moves  past  them.  An  increase  in  the  relative  speed  of  this 
passing  texture  brings  about  an  increase  in  optical  flow,  indicating 
that  an  object  is  approaching.  Conversely,  a  decrease  in  optical  flow 
indicates  that  the  surroundings  are  receding.  The  optical  flow  sensor 
described  here  is  designed  to  allow  MAVs  to  maintain  a  constant  al¬ 
titude  (see  Fig.  8)  and  to  avoid  collision  with  objects  while  in  flight 
(see  Fig.  9). 

To  this  point,  the  altitude  function  of  the  sensor  has  been  the  pri¬ 
mary  focus  of  research.  In  July  of  2000,  a  version  of  the  sensor  was 
placed  on  a  remotely  controlled  model  aircraft  and  flown  by  an  oper¬ 
ator  to  assess  the  ability  of  the  sensor  to  measure  optical  flow  and 
maintain  a  constant  altitude.  During  periods  of  straight  flight  the  air¬ 
craft  was  allowed  to  fly  without  adjustment  by  the  operator.  It  was 
observed  that  the  aircraft  reacted  to  changes  in  proximity  to  the 
ground  and  made  attempts  to  adjust  altitude  to  avoid  crashing  or  fly¬ 
ing  too  high.  The  flight  control  was  crude,  however,  due  to  the  sim 
pie  nature  of  the  aircraft  control  system  employed. 

The  sensor  package  flown  in  the  above  test  was  fabricated  through 
MOSIS<5),  a  prototyping  and  small- volume  production  service  for 
very  large  scale  integration  (VLSI)  circuit  development.  The  pack¬ 
age  included  an  analog  optical  flow  sensor,  produced  on  a  2-2  x  2-2 
mm  silicon  chip  with  a  l-2pm  process;  a  Microchip  Technology 
PIC16C76  programmable  micro-controller;  and  associated  wiring, 
timer,  and  battery.  The  micro-controller  processed  the  analog  output 
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Figure  6(a). 

of  the  flow  sensor  to  calculate  the  optical  flow  rate  and  to  determine 
the  direction  of  motion.  This  configuration  allows  some  flexibility  of 
operation,  because  the  micro-controller  can  be  programmed  to  per¬ 
form  a  range  of  functions  with  the  output  data.  However,  the  added 
weight  and  power  consumption  of  the  micro-controller  is  a  liability. 
The  micro  controller’s  processing  speed  is  also  a  limitation,  though 
new  versions  are  now  available  with  significantly  improved  perfor¬ 
mance. 

The  basic  function  of  the  sensor  is  illustrated  in  Fig.  10.  A1  x  1  ar¬ 
ray  of  photoreceptors  on  the  chip  are  stimulated  by  the  presence  of 
contrasting  light  levels  moving  across  the  visual  field.  Movement  is 
detected  and  recorded  by  a  block  of  feature  detectors,  set  to  look  for 
four  or  six  different  patterns  of  transition  between  the  elements  of 
the  photo-detector  array.  A  set  of  winner  takes  all  (WTA)  circuits  in¬ 
dicate  that  a  transition  has  taken  place.  The  speed  measurement  is 
determined  by  the  time  between  WTA  transitions.  The  micro-con¬ 
troller  controls  the  feature  detector  that  is  to  be  used,  but  cycles 
through  all  six  before  giving  a  measurement  of  the  optical  flow. 
When  all  six  have  been  sampled,  the  average  speed  is  determined. 
This  is  the  speed  that  is  used  to  control  the  aircraft.  The  process  is 
explained  in  more  detail  in  (6). 

As  of  this  writing,  the  only  changes  in  the  package  flown  on  the 
test  in  2000  have  been  modifications  to  correct  errors  made  in 
wiring,  transistor  sizes,  and  connections  on  the  chip  itself.  The  result 
is  that  optical  flow  is  now  measured  as  an  average  of  all  six  feature 
detectors,  not  four  as  in  the  previous  version.  This  should  result  in 
significant  improvement  in  flow  rate  measurement,  and  a  more  con¬ 
sistent  average  value.  The  latest  version  of  the  chip  is  now  correct  in 
its  functional  operation,  as  determined  by  bench  testing.  The  addi¬ 
tion  of  two  extra  feature  detector  inputs  will  require  minor  modifica- 


-5.0  0.0  5.0  10.0  15.0  20.1 

a  (deg.) 

Figure  7(a). 


a  (deg.) 


Figure  7(b). 


0.40  . . 1,11 . I  ■  ■  ■  ■  l 


-5.0  0.0  5.0  10.0  15.0  20.0 

a  (deg.) 


Figure  7(c). 


-20.0  -15.0  -10.0  -5.0  0.0  5.0  10.0  15.0  20.0 

0  (deg.) 


Figure  6(b). 


Figure  7(d). 


KIHIII  .  .  HttHH 


Kellog  etal 


The  nrl  micro  tactical  expendable  (mite)  air  vehicle 


NUMBER 


tions  to  the  code  operating  the  micro-controller. 

Current  research  is  focused  of  developing  a  complete,  on-chip 
sensor  that  computes  the  optical  flow  rate  and  direction  with  hard¬ 
ware  only.  This  eliminates  the  need  for  the  computing  power  of  the 
micro-controller,  reducing  size  and  weight.  Preliminary  designs  have 
been  received  from  MOSIS  and  are  currently  undergoing  bench  test¬ 
ing  and  debugging. 

Three  dimensional  computer  vision  for  navigation  and 
guidance 

Numerous  options  have  been  studied  for  providing  visual  perception 
to  an  MAV.  The  visual  capabilities  that  seem  most  essential  for  au¬ 
tonomous  flight  are  collision  avoidance  and  navigation.  Less  vital, 
but  potentially  enormously  useful,  is  the  efficient  recognition  and  lo¬ 
calisation  of  various  targets.  Novel  technology  has  been  developed 
in  all  of  these  categories.  This  work  consists  of  two  major  parts: 
range-based  vision  for  object  recognition  and  pose  estimation,  and 
monocular  vision  for  navigation  and  collision  avoidance. 

A  wide  variety  of  range-based  imaging  methods,  relying  on  mea¬ 
surement  of  the  distance  between  the  vision  system  and  objects  in  its 
surroundings,  were  considered  for  the  recognition  and  localisation 
task.  A  particularly  suitable  approach  is  multinocular  stereo  using  a 
projected,  pseudo-random  dot  pattern.  In  this  approach,  a  xenon 
flash  tube  is  used  to  project  spots  of  light  onto  the  terrain  around  the 
MAV.  While  the  flashed  dot  pattern  increases  the  detectability  of  the 
vehicle,  this  could  be  mitigated  somewhat  by  the  use  of  an  infrared 
flash.  Two  or  three  cameras  on  the  MAV,  of  known  orientation  and 
separation,  image  these  spots  and  the  results  are  processed  to  obtain 
range  information.  This  system  has  the  properties  of  motion  toler¬ 
ance  and  high  resolution.  All  the  components  can  be  small  and  light. 
The  resulting  range  images,  taken  at  a  few  meters  range,  would  sup¬ 
port  the  use  of  Tripod  Operators^7  and  8),  an  efficient  and  general 
method  for  recognising  and  localising  freeform  surface  shapes  in  6 
DOF.  The  automotive  industry  (Ford  and  Perceptron)  is  testing  Tri¬ 
pod  Operator  technology  for  the  recognition  and  6  DOF  pose  estima¬ 
tion  of  transmission  parts  during  assembly,  with  great  success.  Tri¬ 
pod  Operators  would  provide  the  ability  to  recognise  immediately 
upon  encounter  many  kinds  of  targets.  A  disadvantage  of  range- 
based  vision  methods  is  that  all  of  the  many  kinds  of  triangulation 
imagers  have  low  accuracy  beyond  a  few  meters,  and  present  lidars 
are  too  heavy  for  MAV  use.  Nevertheless,  the  future  application  of 
range-based  technology  is  foreseen  in  the  rapid  and  accurate  recog¬ 
nition  and  localisation  of  known  targets  at  close  range. 

Passive  monocular  vision  has  the  advantage  of  complete  passivity, 
and,  after  equal  design  miniaturisation  efforts,  will  be  lighter  in 
weight  than  a  dot  stereo  system.  Its  development  is  currently  a  first 
priority  because  it  will  allow  the  estimation  of  aircraft  motion  in  up 
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Figure  9. 

to  6  DOF,  thus  supporting  navigation.  Also,  it  can  provide  sparse 
range  for  near  points  to  support  obstacle  avoidance.  NRL’s  approach 
is  to  find  corresponding  features  in  successive  camera  images,  and  to 
deduce  from  these  the  MAV’s  relative  pose,  i.e.  the  relative  position 
in  6  DOF  between  the  aircraft  and  its  surroundings.  Two  methods 
are  under  development,  based  on  horizon  registration  and  point  cor¬ 
respondences,  respectively.  The  first  can  serve  as  the  preprocessor 
for  the  second.  In  all  this  work  computational  efficiency  is  para¬ 
mount,  because  it  is  difficult  to  obtain  a  fast  computer  with  suffi¬ 
ciently  low  weight  and  power.  Experiments  are  progressing  with  a 
15g,  486-level  microprocessor  machine,  running  at  60MHz  at  2-5 
Watts.  A  120  x  160-pixel  camera  is  used  in  conjunction  with  this,  to 
minimise  computation. 

Initial  registration  consists  of  an  efficient  procedure  for  detecting 
the  horizon  edge.  A  specialised  edge  detector  exploits  the  shape  of 
the  light-to-dark  transition  at  the  horizon  in  one  image.  A  single, 
one-dimensional  search  along  the  pixels  of  the  detected  horizon  is 
used  to  find  the  rigid  motion  of  the  image  plane  that  registers  the 
horizon  in  two  consecutive  frames.  This  initial  alignment  allows  the 
estimation  of  the  three  angular  coordinates  of  the  aircraft’  s  frame-to- 
frame  motion.  In  this  application,  frame-to-frame  changes  are  very 
great,  so  that  carefully  designed  matching  procedures  are  necessary 
to  find  corresponding  features  in  two  frames. 

After  the  change  in  the  observed  position  of  the  horizon  has  pro¬ 
vided  an  estimation  of  the  MAV’s  frame-to-frame  change  in  attitude, 
the  subsequent  processes  of  estimating  the  MAV’s  translational  mo¬ 
tion  and  of  refining  the  three  attitude  angles  is  simplified.  Using  a 
carefully  designed  interest  point  detector,  sparse  pointlike  features 
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Figure  1 1 . 


are  detected  in  the  two  successive  images.  These  could  be  the  tips  of 
tree  branches  or  the  corner  of  a  building-features  that  can  be  easily 
and  unambiguously  localised  as  points.  Not  suitable  are  roof  sur¬ 
faces,  road  edges,  tree  limbs,  and  the  like,  or  any  object  with  a  planar 
or  linear  appearance.  With  a  small  number  of  appropriately  scattered 
points  it  is  possible  to  compute  the  ranges  of  all  the  points,  simulta¬ 
neously  with  the  change  in  pose  of  the  aircraft,  although  one  overall 
scale  factor  must  be  deduced,  if  it  is  desired,  from  presumed  speed  or 
from  some  known  scene  structure.  The  analysis  is  done  in  the  con¬ 
text  of  epipolar  geometry  (see  Fig.  11).  A  description  of  this  work  is 
in  Pipitone,  et  aK9).  The  resulting  pose  estimates  are  then  chained  to¬ 
gether  to  maintain  a  cumulative  estimate  of  the  aircraft’s  motion. 
Given  a  known  initial  position  and  attitude  at  launch,  the  MAV  is 
able  to  track  its  progress  analogously  to  odometry  with  a  ground  ro¬ 
bot,  but  in  6  DOF.  The  most  useful  of  the  pose  parameters  appears  to 
be  the  azimuthal  angle,  since  it  provides  the  geographic  bearing  for 
navigation. 


CONTROL  OF  LARGE  GROUPS  OF 
AUTONOMOUS  VEHICLES 

Thus  far,  this  paper  has  not  focused  on  the  problem  of  operating  a 
large  number  of  vehicles  in  conjunction  with  each  other.  In  this  sec¬ 
tion,  scenarios  are  considered  in  which  one  might  need  to  control  the 
behavior  of  swarms  of  vehicles,  also  called  ‘agents.’  For  example, 
multiple  autonomous  MITEs  could  provide  a  distributed  network  of 
roving  sensor  platforms,  maintaining  surveillance  over  a  large  area. 
A  general  framework  is  presented,  which  satisfies  the  objective  of 


Figure  13. 


distributed  control  of  a  large  number  of  agents  that  range  in  scale 
from  neurons,  nanobots,  or  micro-electromechanical  systems 
(MEMS)  to  MAVs  and  satellites.  It  is  assumed  that  an  agent’s  sen¬ 
sors  perceive  the  world  (including  other  agents),  and  an  agent’s  ef¬ 
fectors  make  changes  to  that  agent  or  the  world  (including  other 
agents).  Often,  agents  can  only  sense  and  affect  nearby  agents;  thus 
the  problem  is  usually  one  of  ‘local’  control.  Sometimes  control  is 
also  guided  by  global  constraints  and  interactions. 

One  of  the  biggest  problems  is  that  is  often  unknown  how  to  cre¬ 
ate  the  proper  control  rules  for  the  operating  scenarios.  Agent  behav¬ 
ior  in  these  scenarios  must  satisfy  challenging  constraints.  Not  only 
does  the  desired  global  behavior  of  the  agents  need  to  emerge  from 
local  agent  interactions  (i.e.  self-assembly  or  self-organisation),  but 
there  also  needs  to  be  some  measure  of  fault-tolerance,  i.e.  the  global 
behavior  degrades  very  gradually  if  individual  agents  are  damaged. 
Self-repair  is  also  desirable,  where  a  disrupted  or  damaged  system 
repairs  itself. 

Self-assembly,  fault-tolerance,  and  self-repair  are  precisely  those 
principles  exhibited  by  natural  systems.  Thus,  many  answers  to  the 
problems  of  distributed  control  may  lie  in  the  examination  of  the  nat¬ 
ural  laws  of  physics. 

This  section  proposes  a  general  framework,  called  ‘physicomimet- 
ics,’  for  distributed  control  in  which  physics-inspired  forces  control 
agents.  The  term  ‘physicomimetics’  is  used  because  the  theory  is 
motivated,  but  not  restricted  to,  natural  physical  forces.  Physi¬ 
comimetics  forces  are  not  real  forces;  nevertheless  the  agents  can  act 
as  if  the  forces  are  real.  Thus  the  agent’s  sensors  must  see  enough  to 
allow  it  to  compute  the  virtual  force  to  which  it  is  reacting.  The 
agent’ s  effectors  must  allow  it  to  respond  to  this  perceived  force. 

There  are  several  potential  advantages  to  this  approach.  First,  in 
the  real  physical  world,  collections  of  small  entities  yield  surprising- 


Figure  12. 


Figure  14(a). 


Kellog  et al 


The  nrl  micro  tactical  expendable  (mite  )  air  vehicle 


NUMBER 


ly  complex  behavior  from  very  simple  interactions.  Thus  there  is  a 
precedent  for  believing  that  complex  control  can  be  achieved 
through  simple  local  interactions.  This  is  required  for  very  small 
agents  (such  as  neurons  or  nanobots),  since  their  sensors  and  effec¬ 
tors  will  necessarily  be  primitive.  Also,  since  the  approach  is  largely 
independent  of  the  size  and  number  of  agents,  the  results  should 
scale  well  to  larger  agents  and  larger  sets  of  agents.  Finally,  such 
systems  are  potentially  analysable  by  physics  analysis  methods  (e.g. 
statistical  mechanics). 


Physicomimetics  Framework 

The  motivation  for  this  work  stems  from  a  desire  for  swarms  of 
MAVs  to  perform  various  tasks.  This  approach  treats  agents  as  phys¬ 
ical  particles,  which  could  range  in  size  from  nanobots  to  satellites. 
A  simple  but  realistic  physical  simulation  of  the  particles’  behavior 
was  built.  Particles  exist  in  two  (or  three)  dimensional  worlds  and 
are  considered  to  be  point-masses.  Each  particle  p  has  position  p  = 
( x ,  y)  and  velocity  v  =  (vx ,  v}).  A  discrete- time  approximation  to  the 
continuous  behavior  of  the  particles  is  used,  with  time- step  At.  At 
each  time  step,  the  position  of  each  particle  undergoes  a  perturbation 
A p.  The  perturbation  depends  on  the  current  velocity  Ap  =  v  At.  The 
velocity  of  each  particle  at  each  time  step  also  changes  by  Av.  The 
change  in  velocity  is  controlled  by  the  force  on  the  particle  Av=F  At 
lm ,  where  m  is  the  mass  of  that  particle  and  F  is  the  force  on  that 
particle.  A  frictional  force  is  included,  for  self-stabilisation. 

For  MAVs,  the  initial  conditions  are  similar  to  those  of  a  ‘big 
bang’ -the  MAVs  are  assumed  to  be  released  from  a  canister  dropped 
from  an  aeroplane;  they  then  spread  outwards  and  perform  their  task. 
This  is  simulated  by  using  a  two-dimensional  Gaussian  random  vari¬ 
able  to  initialise  the  positions  of  all  particles  (MAVs).  Velocities  of 
all  particles  are  initialised  to  be  0,  and  masses  are  all  1  (although  the 


framework  does  not  require  this).  An  example  initial  configuration 
for  100  particles  is  shown  in  Fig.  12. 

Given  the  initial  conditions  and  some  desired  global  behavior,  the 
sensors,  effectors,  and  force  F  laws  need  to  be  defined  such  that  the 
desired  behavior  emerges.  This  is  explored  in  the  next  few  sections, 
for  various  tasks.  Note  that  these  are  simulations  only;  the  process 
has  yet  to  be  applied  to  actual  agents. 

Creating  Hexagonal  Lattice  Formations 

The  example  considered  here  is  that  of  a  swarm  of  MAVs  whose 
mission  is  to  form  a  hexagonal  lattice,  which  creates  an  effective 
sensing  grid.  Essentially,  such  a  lattice  will  create  a  virtual  antenna 
or  synthetic  aperture  radar  to  improve  the  resolution  of  radar  images. 
A  virtual  antenna  is  expected  to  be  an  important  future  application  of 
MAVs.  Currently,  the  technology  for  MAV  swarms  (and  swarms  of 
other  micro- vehicles  such  as  micro-satellites)  is  in  the  early  research 
stage.  Nevertheless  we  are  developing  the  control  software  now  so 
that  we  will  be  prepared. 

Since  MAVs  have  simple  sensors  and  primitive  CPUs,  the  goal 
was  to  provide  the  simplest  possible  control  rules  that  require  mini¬ 
mal  sensors  and  effectors.  At  first  blush,  creating  hexagons  would 
appear  to  be  somewhat  complicated,  requiring  sensors  that  can  cal¬ 
culate  range,  the  number  of  neighbors,  their  angles,  etc.  However,  it 
turns  out  that  only  range  information  is  required.  To  understand  this, 
recall  the  old  geometry  lesson  in  which  six  circles  of  radius  R  can  be 
drawn  on  the  perimeter  of  a  central  circle  of  radius  R  (the  fact  that 
this  can  be  done  with  only  a  compass  and  straight-edge  can  be 
proven  with  Galois  theory).  Figure  13  illustrates  this  construction.  If 
the  particles  (shown  as  small  circular  spots)  are  deposited  at  the  in¬ 
tersections  of  the  circles,  they  form  a  hexagon. 

To  map  this  into  a  force  law,  imagine  that  each  particle  repels 
other  particles  that  are  closer  than  R ,  while  attracting  particles  that 
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Figure  14(f). 


are  further  than  R  in  distance.  Thus  each  particle  can  be  considered 
to  have  a  circular  ‘potential  well’  around  itself  at  radius  ^-neighbor¬ 
ing  particles  will  want  to  be  at  distance  R  from  each  other.  The  inter¬ 
section  of  these  potential  wells  is  a  form  of  constructive  interference 
that  creates  ‘nodes’  of  very  low  potential  energy  where  the  particles 
will  be  likely  to  reside  (again  these  are  the  small  circular  spots  in  the 
previous  figure).  Thus  the  particles  serve  to  create  the  very  potential 
energy  surface  to  which  they  respond.  The  entire  potential  energy 
surface  is  never  actually  computed.  Particles  only  compute  local 
force  vectors  for  their  current  locations. 

With  this  in  mind  a  force  law  is  defined 


F-G  mi  rrij  /  r2 


where  F  is  the  magnitude  of  the  force  between  two  particles  i  and  j, 
and  r  is  the  range  between  the  two  particles.  The  ‘gravitational  con¬ 
stant’  G  is  set  at  initialisation.  The  force  is  repulsive  if  r  <  R  and  at¬ 
tractive  if  r>R.  Each  particle  has  one  sensor  that  can  detect  the  range 
to  nearby  particles.  The  only  effector  is  to  be  able  to  move  with  ve¬ 
locity  n.  To  ensure  that  the  force  laws  are  as  local  as  possible  in  na¬ 
ture,  particles  have  a  visual  range  of  only  1-5 R. 

The  initial  universe  of  100  particles  (as  shown  in  Figure  12)  is 
now  allowed  to  run,  using  this  very  simple  force  law.  For  a  radius  R 
of  50  we  have  found  that  a  gravitational  constant  of  G  =  1200  pro¬ 
vides  good  results,  as  can  be  seen  in  Fig.  14. 

The  hexagonal  lattice  is  not  perfect-the  final  perimeter  is  not  a 
hexagon,  although  this  is  not  surprising,  given  the  lack  of  global 
constraints.  However,  many  hexagons  are  clearly  embedded  in  the 
structure  and  the  overall  structure  is  quite  hexagonal.  It  is  also  inter¬ 
esting  to  note  that  each  node  in  the  structure  can  have  multiple  parti¬ 
cles  (i.e.  multiple  particles  can  ‘cluster’  together).  Clustering  was  an 
emergent  property  that  was  not  expected,  and  it  provides  increased 
robust  behavior,  because  the  disappearance  (failure)  of  individual 


Figure  16. 


particles  (agents)  from  a  cluster  will  have  minimal  effect.  This  form 
of  fault- tolerance  is  a  result  of  the  setting  of  G. 

Creating  square  lattices 

Given  the  success  in  creating  hexagonal  lattices,  inspiration  led  to 
the  investigation  of  other  regular  structures.  The  square  lattice  is  an 
obvious  choice,  since  (as  with  hexagons)  squares  will  tile  a  2D 
plane.  The  success  of  the  hexagonal  lattice  hinged  upon  the  fact  that 
nearest  neighbors  are  R  in  distance.  Clearly  this  is  not  true  for 
squares,  since,  if  the  distance  between  particles  along  an  edge  is  R , 
the  distance  along  the  diagonal  is  V2 R.  The  problem  is  that  the  parti¬ 
cles  have  no  way  of  knowing  whether  their  relationship  to  neighbors 
is  along  an  edge  or  along  a  diagonal.  Once  again  it  would  appear  as 
if  the  agents  would  need  to  know  angles  or  the  number  of  neighbors 
to  solve  this  difficulty.  In  fact,  a  much  simpler  approach  will  do  the 
trick.  Suppose  at  creation  that  each  particle  is  given  another  attribute, 
called  ‘spin.’  Half  of  the  particles  are  initialised  to  be  spin  ‘up’, 
while  the  other  half  are  initialised  to  be  spin  ‘down.’  This  ‘coloring’ 
of  the  particles  allows  them  to  form  a  square  lattice.  Particles  of 
spin  are  distance  R  from  each  other,  while  particles  of  like  spin  are 
distance  'llR  from  each  other.  Thus  the  sensors  need  to  be  able  to 
detect  one  more  bit  of  information,  spin.  The  same  force  law  is  used 
as  before.  In  this  case,  however,  the  range  r  is  renormalised  to  be 
rN 2  for  two  particles  having  the  same  spin.  Then  once  again  the 
force  is  repulsive  if  r  <  R  and  attractive  if  r  >  R.  The  only  effector  is 
to  be  able  to  move  with  velocity  v. 

The  initial  universe  of  100  particles  is  allowed  to  run,  using  this 
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Figure  18. 

very  simple  force  law.  The  final  result  is  shown  in  Fig.  15. 

The  physicomimetics  framework  can  easily  be  extended  to  three- 
dimensional  worlds.  For  example,  Fig.  16  shows  a  cubic  lattice. 

Dynamic  behaviours 

In  the  previous  examples  the  MAVs  move  into  formation  and  then 
stop.  Dynamic  behaviours  are  also  of  interest,  where  the  MAVs  con¬ 
stantly  move  through  space.  Two  obvious  tasks  of  interest  are  sur¬ 
veillance  and  perimeter  defense.  For  these  two  tasks  the  MAVs  must 
cover  a  space  while  searching  for  intruders. 

To  achieve  this  the  above  framework  is  used,  but  with  inter-agent 
attraction  turned  off.  The  agents  repel  each  other,  but  are  attracted  to 
intruders  that  they  can  detect.  They  also  require  a  means  of  sensing 
the  limits  of  their  operating  area.  The  net  effect  is  similar  to  a  simu¬ 
lation  of  gas  molecules  in  a  container.  Figures  17  and  18  show 
examples  of  MAVs  performing  surveillance  and  perimeter  defense. 
In  both  figures,  an  intruder  is  depicted  by  a  small  square.  Note  how 
some  of  the  agents  cluster  in  the  square,  thus  ‘attacking’  the  intruder. 
These  snapshots  are  static,  and  it  is  important  to  remember  that  the 
agents  are  in  fact  continually  moving. 

A  nice  aspect  of  this  framework  is  its  robustness.  If  agents  are 
lost,  the  ‘gas’  becomes  more  diffuse,  but  the  agents  still  perform  the 
task,  although  somewhat  less  efficiently.  Moreover,  the  addition  of 
new  agents  is  also  handled  gracefully. 

Future  work  will  focus  on  other  behaviors  and  formations.  As 
applications  arise,  they  will  suggest  new  operating  scenarios  for 
which  physicomimetics  solutions  will  be  derived. 


CURRENT  STATUS  OF  MITE 
DEVELOPMENT 

MITE  vehicles  of  various  sizes  and  weights  have  demonstrated  good 
performance  in  flight.  The  vehicle  is  stable,  yet  is  capable  of  agile 
manoeuvring.  Its  control  and  propulsion  systems  continue  to  be  im 
proved.  Experiments  are  now  underway  to  test  vehicle  flight  control 
through  differential  motor  thrust.  This  may  do  away  with  the  need 
for  movable  flight  control  systems.  Techniques  for  electrical  energy 
storage  are  improving  constantly,  driven  by  the  vast  and  rapidly 
expanding  market  for  hand-held  electronic  devices.  Lithium-polymer 
and  lithium-manganese  dioxide  batteries  and  miniaturised  chemical 
fuel  cells  may  be  practical  for  MAV  use  in  the  near  future. 

As  of  this  writing,  MITE  vehicles  have  only  flown  under  remote 
control  by  a  human  operator.  Prototype  versions  of  the  optic  flow 
sensor  and  the  monocular  motion  stereo  system  have  been  flown  in 


other  testbed  vehicles.  The  optic  flow  sensor  has  demonstrated  both 
collision  avoidance  of  a  single  obstacle,  and  altitude  hold  above 
ground  level  capabilities.  The  sensor  element,  detection  methods, 
and  signal  processing  are  being  refined  to  produce  a  smaller,  more 
reliable  unit.  The  monocular  motion  stereo  system  has  acquired  se¬ 
quential  in-flight  images  that  are  being  used  to  further  develop  the 
image  processing  algorithms.  Also  in  development  is  a  conventional 
autopilot  of  appropriate  size  and  weight  for  the  MITE,  using  com 
mercially  available  sensors.  It  is  envisioned  that  the  vision-based 
sensors  will  be  incorporated  into  this  autopilot  as  they  become  avail¬ 
able. 
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